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Introduction 

Percolation as a mathematical theory was intro- 
duced by Broadbent and Hammersley 0, as a 
stochastic way of modeling the flow of a fluid or 
gas through a porous medium of small channels 
which may or may not let gas or fluid pass. It 
is one of the simplest models exhibiting a phase 
transition, and the occurrence of a critical phe- 
nomenon is central to the appeal of percolation. 
Having truly applied origins, percolation has been 
used to model the fingering and spreading of oil 
in water, to estimate whether one can build non- 
defective integrated circuits, to model the spread 
of infections and forest fires. From a mathemati- 
cal point of view percolation is attractive because 
it exhibits relations between probabilistic and al- 
gebraic / topological properties of graphs. 

To make the mathematical construction of such 
a system of channels, take a graph g (which orig- 
inally was taken as Z''), with vertex set and 
edge set £ , and make all the edges independently 
open (or passable) with probability p or closed (or 
blocked) with probability 1 — p. Write Pp for the 
corresponding probability measure on the set of 
configurations of open and closed edges — that 
model is called bond percolation. The collection of 
open edges thus forms a random subgraph of g , 
and the original question stated by Broadbent was 
whether the connected component of the origin in 
that subgraph is finite or infinite. 

A path on is a sequence v\,V2t ■ ■ of vertices 
of g , such that for all i > 1, c, and are adjacent 
on ^ . A path is called open if all the edges {c, , J^,+i} 
between successive vertices are open. The rnfinite- 
ness of the cluster of the origin is equivalent to 
the existence of an imboimded open path starting 
from the origin. 

There is an analogous model, called site percola- 
tion, in which all edges are assumed being pass- 
able, but the vertices are independently open or 
closed with probability p or 1 — p, respectively. An 
open path is then a path along which all vertices 



are open. Site percolation is more general than 
bond percolation in the sense that the existence of 
a path for bond percolation on a graph g is equiv- 
alent to the existence of a path for site percolation 
on the covering graph of g . However, site perco- 
lation on a given graph may not be equivalent to 
bond percolation on any other graph. 

All graphs imder consideration will be assumed 
to be connected, locally finite and quasi-transitive. 
If A, B C 1^ , then A ^ B means that there exists an 
open path from some vertex of A to some vertex of 
B; by a slight abuse of notation, m <-> c will stand 
for the existence of a path between sites u and v, 
i.e. the event {«} <-> {v}. The open cluster C{v) of 
the vertex v is the set of all open vertices which are 
connected to v by an open path: 



C{v) = {uEl' : 



V}. 



The central quantity of the percolation theory is 
the percolation probability: 

e{p) := Pp{0 ^ oo} = Pp{\C{0)\ = oo}. 



The most important property of the percola- 
tion model is that it exhibits a phase transition, i.e. 
there exists a threshold value pc G [0, 1], such that 
the global behavior of the system is substantially 
different in the two regions p < pc and p > pc. 
To make this precise, observe that 6* is a non- 
decreasing function. This can be seen using Ham- 
mersley's joint construction of percolation systems 
for all p e [0, 1] on ^: Let {U{v), e 1^ } be inde- 
pendent random variables, uniform in [0,1]. De- 
clare V to be p-open if U{v) < p, otherwise it is 
declared p-closed. The configuration of p-open 
vertices has the distribution Pp for each p G [0, 1]. 
The collection of p-open vertices is non-decreasing 
in p, and therefore 6{p) is non-decreasing as well. 
Clearly 61(0) = and 61(1) = 1. 

The critical probability is defined as 

Pc ■■= Pc{g) = sup{p : e{p) = 0}. 
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Figure 1: The behavior of 6'(p) around the critical 
point (for bond-percolation) 

By definition, when p < Pc the open cluster of the 
origin is Pp-a.s. finite, hence all the clusters are also 
finite. On the other hand, for p > pc there is a 
strictly positive Pp-probability that the cluster of 
the origin is infinite. Thus, from Kolmogorov's 
zero-one law it follows that 

Pp{|C(u)| = oo for some v E 1^} = 1 for p > pc- 

Therefore, if the intervals [0, pc) and {pc, 1] are both 
non-empty, there is a phase transition at pc- 

Using a so-called Peierls argument it is easy to 
see that pc{(^) > for any graph g of boimded de- 
gree. On the other hand, Hammersley proved that 
pd^'') < 1 for bond percolation as soon as d >2, 
and a similar argument works for site percolation 
and various periodic graphs as well. But for some 
graphs g it is not so easy to show that pdg) < 1. 
One says that the system is in the subcritical (resp. 
supercritical) phase if p < (resp. p > pc). 

It was one of the most remarkable moments in 
the history of percolation when Kesten proved flOf 
that the critical parameter for bond-percolation on 
1? is equal to 1/2. Nevertheless the exact value of 
pdg) is known only for a handful of graphs, all of 
them periodic and two-dimensional — see below. 

1 Percolation in l/ 

The graph on which most of the theory was origi- 
nally built is the cubic lattice Z'', and it was not be- 
fore the late 20th century that percolation was se- 
riously considered on other kinds of graphs (such 
as e.g. Cayley graphs), on which specific phenom- 
ena can appear, such as the coexistence of multi- 
ple infinite clusters for some values of the param- 
eter p. In all this section, the underlying graph is 
thus assumed to be l/ for d >2, although most 
of the results still hold in the case of a periodic d- 
dimensional lattice. 



1.1 The sub-critical regime 

When p < Pc, all open clusters are finite almost 
surely. One of the greatest challenges in per- 
colation theory has been to prove that xip) ■= 
£p{|C(c)|} is finite \t p < pc {Ep stands for the ex- 
pectation with respect to Pp). For that one can 
define another critical probability as the threshold 
value for the finiteness of the expected cluster size 
of a fixed vertex: 

Prig) :=sup{p: xip) < oo}. 

It was an important step in the development of 
the theory to show that prig) = pd(^)- The fun- 
damental estimate in the subcritical regime, which 
is a much stronger statement than prig) = pdQ)' 
is the following: 

Theorem 1 (Aizenman and Barsky, Menshikov) 

Assume that g is periodic. Then for p < Pc there exist 
constants < Ci , C2 < oo, such that 

Pp{\C{v)\ > n} < Cie-^^" . 

The last statement can be sharpened to a "local 
limit theorem" with the help of a subadditivity ar- 
gument : For each p < p^ there exists a constant 
< Csip) < 00, such that 

Jim-llogPp{|C(i;)|=n} = C3(p). 

1.2 The super-critical regime 

Once an infinite open cluster exists, it is natural 
to ask how it looks like, and how many infinite 
open clusters exist. It was shown by Newman and 
Schulman that for periodic graphs, for each p, ex- 
actly one of the following three situations prevails: 
If N G Z+ U {00} is the number of infinite open 
clusters, then Pp(N = 0) = 1, or P/N = 1) = 1, or 
Pp(N = 00) = 1. 

Aizenman, Kesten and Newman showed that 
the third case is impossible on Z"^. By now several 
proofs exist, perhaps the most elegant proof of that 
is due to Burton and Keane, who prove that indeed 
there cannot be infinitely many infinite open clus- 
ters on any amenable graph. However, there are 
some graphs, such as regular trees, on which coex- 
istence of several infinite clusters is possible. 

The geometry of the infinite open cluster can be 
explored in some depth by studying the behavior 
of a random walk on it. When d = 2, the random 
walk is recurrent, and when d > 3 is a.s. transient. 
In all dimensions d > 2 the walk behaves diffu- 
sively, and the Central Limit Theorem and the In- 
variance principle were established in both the an- 
nealed and quenched cases. 
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Wulff droplets 

In the supercritical regime, aside from the infinite 
open cluster, the configuration contains finite clus- 
ters of arbitrary large sizes. These large finite open 
clusters can be thought of as droplets swimming 
in the areas surrounded by an infinite open cluster. 
The presence at a particular location of a large fi- 
nite cluster is an event of low probability, namely, 
on l/, d > 2, for p > pc, there exist positive con- 
stants < Ci{p), Csip) < oo, such that 



Ciip) < 



1 



„(d-i)/d 



logPp{|C(z;)| = «}<C5(p) 



for all large n. This estimate is based on the fact 
that the occurrence of a large finite cluster is due to 
a surface effect. The typical structure of the large 
finite cluster is described by the following theo- 
rem: 

Theorem 2 Let d > 2, and p > pc- There exists a 
hounded, closed, convex subset W of R** containing 
the origin, called the normalized Wulff crystal of the 
Bernoulli percolation model, such that, under the condi- 
tional probability Pp{- I n'^ < |C(0)| < oo} , the random 
measure ^ 

" .TeC(o) 

(where 5x denotes a Dirac mass at x) converges weakly 
in probability towards the random measure ^(p)lw(x — 
M)dx (where M is the rescaled center of mass of the 
cluster C(0)j. The deviation probabilities behave as 
exp{— cn'*"^} (i.e. they exhibit large deviations of sur- 
face order). 

This result was proved in dimension 2 by 
Alexander, Chayes and Chayes yj, and in dimen- 
sions 3 and more by Cerf 

1.3 Percolation near the critical point 
1.3.1 Percolation in slabs 

The main macroscopic observable in percolation is 
0{p), which is positive above pc, below pc, and 
continuous on [0, 1] \ {pc}- Continuity at pc is an 
open question in the general case; it is known to 
hold in two dimensions (cf. below) and in high 
enough dimension (at the moment d > 19 though 
the value of the critical dimension is believed to be 
6) using lace expansion methods. The conjecture 
that 9{pc) = for 3 < d < 18 remains one of the 
major open problems. 

Efforts to prove that led to some interesting and 
important results. Barsky, Grimmett and New- 
man solved the question in the half-space case. 



and simultaneously showed that the slab perco- 
lation and half-space percolation thresholds coin- 
cide. This was complemented by Grimmett and 
Marstrand showing that 

p,{slab) = Pc{Z'). 
1.3.2 Critical exponents 

In the sub-critical regime, exponential decay of the 
correlation indicates that there is a finite correlation 
length ^(p) associated to the system, and defined 
(up to constants) by the relation 



Pp(0 nx) w exp 



nip{x) 



where (p is boimded on the unit sphere (this is 
known as Ornstein-Zernike decay). The phase tran- 
sition can then also be defined in terms of the di- 
vergence of the correlation length, leading again 
to the same value for p^-; the behavior at or near 
the critical point then has no finite characteristic 
length, and gives rise to scaling exponents (conjec- 
turally in most cases). 

The most usual critical exponents are defined 
as follows, if 9{p) is the percolation probability, C 
the cluster of the origin, and ^(p) the correlation 
length: 
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These exponents are all expected to be univer- 
sal, i.e. to depend only on the dimension of the 
lattice, although this is not well understood at the 
mathematical level; the following scaling relations 
between the exponents are believed to hold: 

2-a = 7 + 2/3 = /3(5 + l), A = 5/3, 7 = ^2-?]). 

In addition, in dimensions up to dc = 6, two addi- 
tional hyperscaling relations involving d are strongly 
conjectured to hold: 



dp = 5 + \, dv = 2 



a. 
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while above dc the exponents are believed to take 
their mean-field value, i.e. the ones they have for 
percolation on a regular tree: 

a = -l, (3 = 1, 7 = 1, 5 = 2, 

1 1 

r? = 0, 1/ = p= -, A = 2. 

Not much is known rigorously on critical ex- 
ponents in the general case. Hara and Slade 
(i9i) proved that mean field behavior does hap- 
pen above dimension 19, and the proof can likely 
be extended to treat the case d > 7. In the two- 
dimensional case on the other hand, Kesten ( llSl ) 
showed that, assuming that the exponents 5 and 
p exist, then so do /3, 7, rj and v, and they satisfy 
the scaling and hjrperscaling relations where they 
appear. 

1.3.3 The incipient infinite cluster 

When studying long-range properties of a critical 
model, it is useful to have an object which is infi- 
nite at criticality, and such is not the case for per- 
colation clusters. There are two ways to condition 
the cluster of the origin to be infinite when p = pc'. 
The first one is to condition it to have diameter at 
least n (which happens with positive probability) 
and take a limit in distribution as n goes to infin- 
ity; the second one is to consider the model for 
parameter p > pc, condition the cluster of to be 
infinite (which happens with positive probability) 
and take a limit in distribution as p goes to pc. The 
limit is the same in both cases, it is known as the 
incipient infinite cluster. 

As in the super-critical regime, the structure of 
the cluster can be investigated by studying the be- 
havior of a random walk on it, as was suggested by 
de Gennes; Kesten proved that in two dimensions, 
the random walk on the incipient infinite cluster 
is sub-diffusive, i.e. the mean square displacement 
after n steps behaves as «^^^ for some £ > 0. 

The construction of the incipient infinite clus- 
ter was done by Kesten in two dimensions |12j, 
and a similar construction was performed recently 
in high dimension by Van der Hofstad and Jarai 

(iQi)- 

2 Percolation in two dimensions 

As is the case for several other models of statistical 
physics, percolation exhibits many specific prop- 
erties when considered on a two-dimensional lat- 
tice: Duality arguments allow for the computation 



of Pc in some cases, and for the derivation of a pri- 
ori bounds for the probability of crossing events 
at or near the critical point, leading to the fact that 
G{Pc) = 0. On another front, the scaling limit of crit- 
ical site-percolation on the two-dimensional trian- 
gular lattice can be described in terms of SLE pro- 
cesses. 




Figure 2: Two large critical percolation clusters in 
a box of the square lattice (first: bond-percolation, 
second: site-percolation) 



2.1 Duality, exact computations and 
RSW theory 

Given a planar lattice L, define two associated 
graphs as follows. The dual lattice l' has one ver- 
tex for each face of the original lattice, and an edge 
between two vertices if and only if the correspond- 
ing faces of L share an edge. The star graph L* is 
obtained by adding to L an edge between any two 
vertices belonging to the same face (x * is not pla- 
nar in general; {£,£*) is commonly known as a 
matching pair). Then, a result of Kesten is that, un- 
der suitable technical conditions, 

pbond(^) + p^°"'l(x') = pf'iL) + pf" (X*) = 1. 

Two cases are of particular importance: The lat- 
tice is isomorphic to its dual; the triangular lat- 
tice 'T is its own star graph. It follows that 

p^nZ^) = pf^(T) = l. 
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The only other critical parameters that are known 
exactly are pj?°"'*('r) = 2sin(7r/18) (and hence also 
pbond fgj. ^ I ^ g hexagonal lattice), and ff^'^ 
for the bow-tie lattice which is a root of the equa- 
tion — 6p^ + 6p^ + p — 1 = 0. The value of the 
critical parameter for site-percolation on might 
on the other hand never be known, it is even pos- 
sible that it is "just a number" without any other 
signification. 

Still using duality, one can prove that the prob- 
ability, for bond-percolation on the square lattice 
with parameter p = 1/2, that there is a connected 
component crossing an (n + 1) x n rectangle in the 
longer direction is exactly equal to 1/2. This and 
clever arguments involving the symmetry of the 
lattice lead to the following result, proved inde- 
pendently by Russo and by Seymour and Welsh 
and known as the RSW theorem: 

Theorem 3 (Russo 1 15J; Seymour-Welsh ^) For 

every a,b > there exist r] > and Mq > suc/z that 
for every n > Hq, the probability that there is a cluster 
crossing an [na\ x [nb\ rectangle in the first direction 
is greater than rj. 

The most direct consequence of this estimate is 
that the probability that there is a cluster going 
aroimd an annulus of a given modulus is boimded 
below independently of the size of the annulus; 
in particular, almost surely there is some annulus 
aroimd in which this happens, and that is what 
allows to prove that 0{pc) = for bond-percolation 
on I?. 

2.2 The scaling limit 

RSW-type estimates give positive evidence that a 
scaling limit of the model should exist; it is in- 
deed essentially sufficient to show convergence of 
the crossing probabilities to a non-trivial limit as n 
goes to infinity The limit, which should depend 
only on the ration a/b, was predicted by Cardy us- 
ing conformal fields theory methods. A most cel- 
ebrated result of Smirnov is the proof of Cardy's 
formula in the case of site-percolation on the trian- 
gular lattice T : 

Theorem 4 (Smirnov 1181 ) Let Q. be a simply con- 
nected domain of the plane with four points a, b, c, d (in 
that order) marked on its boundary. For every 5 > 
consider a critical site-percolation model on the inter- 
section ofQ. with 5l and let fs{ab, cd;Q.) be the prob- 
ability that it contains a cluster connecting the arcs ab 
and cd. Then: 

1. fs{ab, cd; £2) has a limit fo{ab, cd; Q.)as 5 ^ 0; 



2. The limit is conformally invariant, in the follow- 
ing sense: If ^ is a conformal map from £1 to 
some other domain Q.' = and maps a to a' , 
b to b', c to c' and d to d' , then fo{ab,cd;Q.) = 
fo{a'b',c'd';Q.'); 

3. In the partictdar case when £2 zs an equilateral tri- 
angle of side length 1 and vertices a, b and c, and 
if d is on (ca) at distance x G (0, \) from c, then 
fo{ab, cd;Q.) = x. 

Point 3. in particular is essential since it allows to 
compute the limiting crossing probabilities in any 
conformal rectangle. In the original work of Cardy, 
he made his prediction in the case of a rectangle, 
for which the limit involves hypergeometric func- 
tions; the remark that the equilateral triangle gives 
rise to nicer formulae is originally due to Carleson. 

To precisely state the convergence of percola- 
tion to its scaling limit, define the random curve 
known as the percolation exploration path (see fig.|3) 
as follows: In the upper half -plane, consider a site- 
percolation model on a portion of the triangular 
lattice and impose the boundary conditions that 
on the negative real half -line all the sites are open, 
while on the other half-line the sites are closed. 
The exploration curve is then the common bound- 
ary of the open cluster spanning from the nega- 
tive half -line, and the closed cluster spanning from 
the positive half-line; it is an infinite, self-avoiding 
random curve in the upper half-plane. 




Figure 3: A percolation exploration path 



As the mesh of the lattice goes to 0, the ex- 
ploration curve then converges in distribution to 
the trace of an SLE process, as introduced by 
Schramm, with parameter k = 6 — see fig. |4] The 
limiting curve is not simple anymore (which cor- 
responds to the existence of pivotal sites on large 
critical percolation clusters), and it has Hausdorff 
dimension 7/4. For more details on SLE processes, 
see e.g. the related entry in the present volume. 
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3.1 Percolation 
graphs 



on non-amenable 



Figure 4: An SLE process with parameter k = 6 
(infinite time, with the driving process stopped at 
time 1) 



As an application of this convergence result, one 
can prove that the critical exponents described in 
the previous section do exist (still in the case of the 
triangular lattice), and compute their exact values, 
except for a, which is still listed here for complete- 
ness: 



a 



2 

"3' 
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43 ._91 

18' y 

48 , 91 

y'^=36- 



These exponents are expected to be universal, in the 
sense that they should be the same for percolation 
on any two-dimensional lattice; but at the time of 
this writing this phenomenon is far from being un- 
derstood on a mathematical level. 

The rigorous derivation of the critical exponents 
for percolation is due to Smirnov and Werner 1191 : 
the dimension of the limiting curve was obtained 
by BeffaraEI. 

3 Other lattices and percolative 
systems 

Some modifications or generalizations of standard 
Bernoulli percolation on Z'' exhibit an interesting 
behavior and as such provide some insight into the 
original process as well; there are too many mathe- 
matical objects which can be argued to be percola- 
tive in some sense to give a full account of all of 
them, so the following list is somewhat arbitrary 
and by no means complete. 



The first modification of the model one can think 
of is to modify the imderlying graph and move 
away from the cubic lattice; phase transition still 
occurs, and the main difference is the possibil- 
ity for infinitely many infinite clusters to coexist. 
On a regular tree, such is the case whenever p G 
{pc 1), the first non-trivial example was produced 
by Grimmett and Newman as the product of Z by 
a tree: There, for some values of p the infinite clus- 
ter is unique, while for others there is coexistence 
of infinitely many of them. The corresponding def- 
inition, due to Benjamini and Schramm, is then the 
following: If N is as above the number of infinite 
open clusters, 

p„:=inf{p:P,,(N = l) = l} >p<:. 

The main question is then to characterize graphs 
on which < pc < p,, < 1. 

A wide class of interesting graphs is that of Cay- 
ley graphs of infinite, finitely generated groups. 
There, by a simultaneous result by Haggstrom and 
Peres and by Schonmann, for every p G {pcPu) 
there are Pp-a.s. infinitely many infinite cluster, 
while for every p G (pu , 1] there is only one — note 
that this does not follow from the definition since 
new infinite components could appear when p is 
increased. It is conjectured that pc < Pu for any 
Cayley graph of a non-amenable group (and more 
generally for any quasi-transitive graph with pos- 
itive Cheeger constant), and a result by Pak and 
Smirnova is that every infinite, finitely generated, 
non-amenable group has a Cayley graph on which 
Pc > Pit', this is then expected not to depend on the 
choice of generators. In the general case, it was 
recently proved by Gaboriau that if the graph g 
is unimodular, transitive, locally finite and sup- 
ports non-constant harmonic Dirichlet functions 
{i.e. harmonic fimctions whose gradient is in l^), 
then indeed pdQ) < Pu{§)- 

For reference and further reading on the topic, 
the reader is advised to refer to the review paper 
by Benjamini and Schramm 13J, the lecture notes of 
Peres IT4l . and the more recent article of Gaboriau 

m 

3.2 Gradient percolation 

Another possible modification of the original 
model is to allow the parameter p to depend on 
the location; the porous medium may for instance 
have been created by some kind of erosion, so that 
there will be more open edges on one side of a 
given domain than on the other. If p still varies 
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smoothly, then one expects some regions to look 
subcritical and others to look supercritical, with 
interesting behavior in the vicinity of the critical 
level set {p = pc}. This particular model was in- 
troduced by Sapoval et al. 1161 under the name of 
gradient percolation; see fig. 1^1 




Figure 5: Gradient percolation in a square. In black 
is the cluster spanning from the bottom side of the 
square. 

The control of the model away from the critical 
zone is essentially the same as for usual Bernoulli 
percolation, the main question being how to esti- 
mate the width of the phase transition. The main 
idea is then the same as in scaling theory: If the 
distance between a point v and the critical level 
set is less than the correlation length for parameter 
p{v), then V is in the phase transition domain. This 
of course makes sense only asymptotically, say in 
a large n x n square with p{x, y) = 1 — y/n as is the 
case in the figure: The transition then is expected 
to have width of order n" for some exponent a > 0. 



3.3 First-passage percolation 

First passage percolation (also known as Eden or 
Richardson model) was introduced by Hammer- 
sley and Welsh in 1965 as a time-dependent 
model for the passage of fluid through a porous 
medium. To define the model, with each edge e 6 
E (Z"*) is associated a random variable T(e), which 
can be interpreted as being the time required for 
fluid to flow along e. T{e) are assumed to be in- 
dependent non-negative random variables having 
common distribution F. For any path tt we define 



the passage time r(7r) of tt as 

T(7r) := ^ Tie). 

The first passage time a{x, y) between vertices x and 
y is given by 

y) = rnf {r(7r) : tt a path from x to y}; 

and we can define 

W(0 := {xeZ'':fl(0,x)<0, 

the set of vertices reached by the liquid by time t. 
It turns out that W(f) grows approximately linearly 
as time passes, and that there exists a non-random 
limit set B such that either B is compact and 

(1 - e)B C 1 W(f) C (1 + e)B, eventually a.s. 

for all e > 0, or B = K"*, and 

1 ~ 

{xeM.'' : \x\ < X} C -]N(t), eventually a.s. 

for all X > 0. Here = {2 + {-1/1,1 /if : z e 
W(f)}. 

Studies of first passage percolation brought 
many fascinating discoveries, including King- 
man's celebrated sub-additive ergodic theorem. In 
recent years interest has been focused on study of 
fluctuations of the set W(f) for large t. In spite of 
huge effort and some partial results achieved, it 
still remains a major task to establish rigorously 
conjectures predicted by Kardar-Parisi-Zhang the- 
ory about shape fluctuations in first passage per- 
colation. 

3.4 Contact processes 

Introduced by Harris and conceived with biolog- 
ical interpretation, the contact process on Z'' is a 
continuous-time process taking values in the space 
of subsets of Z''. It is informally described as fol- 
lows: Particles are distributed in Z'' in such a way 
that each site is either empty or occupied by one 
particle. The evolution is Markovian: Each parti- 
cle disappears after an exponential time of param- 
eter 1, independently from the others; at any time, 
each particle has a possibility to create a new parti- 
cle at any of its empty neighboring sites, and does 
so with rate A > 0, independently of everything 
else. 

The question is then whether, starting from a fi- 
nite population, the process will die out in finite 
time or whether it will survive forever with posi- 
tive probability. The outcome will depend on the 
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value of A, and there is a critical value \c, such that 
for A < Ac process dies out, while for A > Ac indeed 
there is survival, and in this case the shape of the 
population obey a shape theorem similar to that of 
first-passage percolation. 

The analogy with percolation is strong, the cor- 
responding percolative picture being the follow- 
me: In Z'l+\ each edge is open with probability 
p e (0, 1), and the question is whether there exists 
an infinite oriented path vr {i.e. a path along which 
the sum of the coordinates is increasing), com- 
posed of open edges. Once again, there is a critical 
parameter customarily denoted by pc, at which no 
such path exists (compare this to the open ques- 
tion of the continuity of the function 6 at pc in di- 
mensions 3 < d < 18). This variation of percolation 
lies in a different universality class than the usual 
Bernoulli model. 

3.5 Invasion percolation 

Let X(e) : e e £ be independent random variables 
indexed by the edge set £ of Z'', d >2, each hav- 
ing uniform distribution in [0,1]. One constructs 
a sequence C = {Q, / > 1} of random connected 
subgraphs of the lattice in the following iterative 
way: The graph Q contains only the origin. Hav- 
ing defined C„ one obtains C,+i by adding to C, an 
edge e,+i (with its outer lying end-vertex), chosen 
from the outer edge boimdary of C, so as to min- 
imize X(e,_|_i). Still very little is known about the 
behavior of this process. 

An interesting observation, relating 0{pc) of 
usual percolation with the invasion dynamics, 
comes from CM. Newman: 

d{pc) = <^ P{x e C} ^ as |x| 00. 

Further reading 

For a much more in-depth review of percolation on 
lattices and the mathematical methods involved in 
its study, and for the proofs of most of the results 
we could only point at, we refer the reader to the 
standard book of Grimmett [7]; another excellent 
general reference, and the only place to find some 
of the technical graph-theoretical details involved, 
is the book of Kesten flT^. More information in the 
case of graphs that are not lattices can be f oimd in 
the lecture notes of Peres 1141 . 

For curiosity, the reader can refer to the first 
mention of a problem close to percolation, in the 
problem section of the first volume of the Amer- 
ican Mathematical Monthly ||2T] . References on 



more specific topics are given at the end of each 
section. 

See also 

Introductory article: Statistical mechanics; 2D 
Ising model; Wulff droplets; Stochastic Loewner 
evolutions. 
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